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ABSTRACT 

We investigate where brightest cluster galaxies (BCGs) sit on the fundamental plane 
of black hole (BH) activity, an established relation between the X-ray luminosity, the 
radio luminosity and the mass of a BH. Our sample mostly consists of BCGs that 
lie at the centres of massive, strong cooling flow clusters, therefore requiring extreme 
mechanical feedback from their central active galactic nucleus (AGN) to offset cooling 
of the intracluster plasma (-L mec h > 10 44-45 erg s _1 ). Based on the BH masses derived 
from the A/bh — cr and A/bh — Mk correlations, we find that all of our objects arc offset 
from the plane such that they appear to be less massive than predicted from their X- 
ray and radio luminosities (to more than a 99 per cent confidence level). For these 
objects to be consistent with the fundamental plane, the Mbh — cr and A/bh — Mk 
correlations therefore seem to underestimate the BH masses of BCGs, on average by 
a factor of 10. Our results suggest that the standard relationships between BH mass 
and host galaxy properties no longer hold for these extreme galaxies. Furthermore, our 
results imply that if these BHs follow the fundamental plane, then many of those that 
lie in massive, strong cool core clusters must be ultramassive with Mbh > 10 10 M Q . 
This rivals the largest BH masses known and has important ramifications for our 
understanding of the formation and evolution of BHs. 
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1 INTRODUCTION 

Brightest cluster galaxies (BCGs) are the most massive 
galaxies in the present-day Universe. They lie at the very 
centres of galaxy clusters and exhibit some of the rich- 
est phenomena known. Many have an active galactic nu- 
cleus (AGN) at their core that is capable of inflating large 
cavities of relativistic plasma through jetted outflows. The 
energetics of these outflows often exceed 10 44-45 erg s _1 
in massive clusters and are sufficient, in the majority of 
clusters, to prevent catast rophic cooling of the surround- 
ing hot X-ray emitting gas dFabian et al|2003l : iBirzan et alJ 
20041 : iForman et al.ll2005l: iRafferty et al-lhood; iFabian et all 
2009: iDunn & Fabian 20061. 12008 : ISanders fc Fabianl I2007F 
McNamara fc N ulscn 20(53). 

Although much is known about the properties of the 
outflows, the fundamental properties of the black holes 
(BHs) powering these outbursts remain largely unknown. 
Since BCGs are the most massive (and luminous) galaxies in 
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the local Universe, by simple extrapolation of observed cor- 
relations, they should host the most massive BHs. Yet, only 
a handful have reliable BH mass measurements from dynam- 
ical modelling of the kinematics: M87 in Virgo (Macchetto 
et al. 1997), NGC1399 in Fornax (Houghton et al. 2006), 
NGC 6086 in A2162 (McConnell et al. 2011b), the BCGs in 
A3565, A1836 and A2052 (Dalla Bonta et al. 2009), as well 
as NGC 3842 in A1367 and NGC 4889 in Coma (McConnell 
et al. 2011a). 

Furthermore, BCGs are galaxies that lie in some of the 
most extreme environments, subject to major mergers in 
the past, and powerful AGN feedback at present times. It is 
therefore not clear if they follow the standard Mbh — cr or 
Mbh — Mk correlations observed in the lower mass galaxies 
(see von der Linden et al. 2007, Dalla Bonta et al. 2009, and 
especially Lauer et al. 2007). The few with reliable mass es- 
timates suggest that some may follow the correlations, while 
others are offset such that the BH mass measured from dy- 
namical modelling is larger than the values predicted from 
the Mbh - cr and A/ B h - M K correlations (e.g. M87, A1836- 
BCG, NGC 3842 and NGC 4889). McConnell et al. (2011a) 
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have also found the most massive BH to date (NGC 4889 
with a 2.1 x 1O 1O M0 BH), providing the first direct evidence 
for the existence of ultramassive BHs (UMBHs, hereafter 
BHs with masses exceeding 10 10 Mq). 

The existence of UMBHs in BCGs has however already 
been predicted, especially for BCGs that lie in the most 
massive and strong cool core clusters (Lx > 10 45 erg s _1 ; 
tcooi < 3 Gyr). In these clusters, the central BH must be 
providing extreme amounts of energy to prevent the sur- 
rounding gas from cooling such that L mcc h > 10 45 erg s~ 4 . 
At these levels, the BH power excee ds 1 per cent of the Ed - 
dington luminosity for a 1O 9 M BH. IChurazov et al.1 (120051 ) 
have argued, by analogy with the X-ray binaries, that BHs 
operating at such powers must be radiatively efficient, while 
at low powers (less than 1 per cent of the Eddington lu- 
minosity), their radiative efficiency drops steadily and the 
power is increasingly taken up by outflows. This explains 
the behaviour of BHs in low/medium mass galaxy clus- 
ters, but cannot explain those in the most massive and 
strong cool core clusters with > 10 45 erg/s jet powers. Since 
they are operating at powers exceeding 1 per cent of the 
Eddington luminosity, they should be radiatively efficient 
(i.e. we should see an X-ra y poin t source). However, in 
iHlavacek-Larrondo &i Fabianl (|201lf l we showed that many 
of these extreme clusters showed no evidence of an X-ray 
point source at their centres. For them to be consistent with 
the lower mass clusters, we proposed that the BCGs host an 
UMBH at their centres (> 10 10 Mq). In this case, the power 
only exceeds 0.1 per cent of the Eddington luminosity and 
they do not require radiatively efficient nuclei. The existence 
of an UMBH has also been proposed to explain the unusu- 
ally large pair of AGN-driven outflows in MS0735.6+7421 
that r equire extreme jet pow ers of 10 46 erg s _1 to create 
them (|McNamara et alj[2009h . 



Here, 



extend 



the 



work 



of 



IHlavacek-Larrondo fc Fabianl |201ll ) and investigate 
where BCGs sit on the fundamental plane (FP) of BH 
activity. Our results show that the A/bh — o and A/bh — Mk 
correlations for BCGs systematically underestimate the 
BH masses, and that many of the BCGs lying in massive 
and strong cool core clusters must have an UMBH at their 
centres with Mbh > 10 10 Mq. In Section 2, we present 
the sample and observations used to determine where the 
objects sit in the FP, and then in Section 3, we present 
the results. Finally, in Section 4 we discuss the results, and 
conclude in Section 5. We adopt Ho = 70kms _1 Mpc 
with f2 m = 0.3, = 0.7 throughout this paper. All errors 
are 2a unless otherwise noted. 



2 FUNDAMENTAL PLANE OF BLACK HOLE 
ACTIVITY 

2.1 Previous studies 

The FP of BH activity is an established correlation relating 
the mass of a compact object to its 2 — 10 keV intrinsic X-ray 
lu minosity (Lx) and 5 GHz core radio lum i nosity (Lsghz), 



ilMerloni et al.1 j2003MFalcke et al.1 (120041) ;|Kording et all 
(|2006l ); iGiiltekin et al.1 (|2009l ); iPlotkin et all (|2012l V (here- 
after MHD2003, FKM2004, KFC2006, GCM2009 and 
PMK2012 respectively). The FP covers over six orders of 



magnitude in mass and is especially important since it pro- 
vides a unification scheme for all BHs and X-ray binaries 
(XRB), regardless of their mass. 

MHD2003 compiled the first detailed study of the FP 
based on a large sample comprising around 100 AGNs and 
8 galactic BHs. The sample consisted of a diverse pop- 
ulation of objects, including X-ray Binaries (XRB), low- 
luminosity AGN (LLAGN), low ionization nuclear emission 
region (LINER), type 1 and 2 Seyfert galaxies, Fanaroff- 
Riley (FR) radio galaxies, as well as radio loud and radio 
quiet quasars. Only direct measurements of Lsghz and Lx 
were used. This study therefore finds the existing correlation 
of the FP averaged over all types of BHs, regardless of the 
accretion state. The best-fitting relationship found by the 
authors is given in Eq. [T] 

logisGHz = 0.601ogLx + 0.781ogM BH + 7.33 (1) 

On the other hand, FKM2004 argued that the FP is 
due to synchrotron emission arising from a jetted BH, and 
therefore only BHs in the "low/hard" (LH) state should be 
included i n the analysi s since this state is dominated by jet 
emission (|Fenderl|200ll ). More recently, both the studies by 
MHD2003 and FKM2004 have been revised by KFC2006, 
with the aim of improving the parameter estimates of the 
FP. Their study suggested that the resulting parameters de- 
pend strongly on the weights given to each AGN class, and 
on the assumptions made on the sources of scatter (e.g. rela- 
tivistic beaming and non-simultaneous measurements of X- 
ray/radio luminosities). They also found that sub-Eddington 
sources equivalent to the LH state seem to follow the plane 
more tightly. In an effort to refine the scatter in the FP, 
GCM2009 only considered sources where dynamical mass- 
measurements were available (18 objects in total). Although 
this limited the sample to nearby AGN (< 30 Mpc), the au- 
thors found that when the sample was further limited to 
low accreting sources (from logA = log(Lx/LBdd) = —4.2 
to logA = —5.2), the scatter in the best-fitting relation de- 
creased, suggesting once more that high accreting sources 
may not belong to the FP. The majority of BCGs we con- 
sider here, are those that lie in massive, strong cool core 
clusters with powerful outflows (L mec h > 10 44- 45 erg s _1 ) 
and radiatively inefficient nuclei (Lx ~ 10 42 erg s _1 ). The 
majority are therefore highly sub-Eddingtion (logA ^ —5, 
see Table [T]), and are considered to be in a state equivalent 
to the LH state. 

PMK2012 have used a more sophisticated regression 
technique (Bayesian), allowing them to further constrain the 
FP and obtain smaller uncertainties on their best-fitting co- 
efficients. They use a sample of low-accreting BHs and find 
that their best-fitting parameters favour the coefficients pre- 
dicted for X-rays that are dominated by optically thin jet 
emission, and not coronal emission, if their objects have 
flat/inverted radio emission. However, they also find that 
the most massive BHs (^10 8 Mq) could be strongly affected 
by synchrotron cooling (see also FKM2004 and KFC2006). 
We discuss this issue further in Section 3.2. 

Our aim is to determine where BCGs sit on the FP 
of BH activity for the general population of compact ob- 
jects. We therefore initially adopt the original correlation 
from MHD2003 to illustrate our results, since it includes all 
sources regardless of the state in which they are in. However, 
the best-fitting relations found by the more recent studies 
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Table 1. Sample of Clusters - (1) Name; (2) redshift; (3) Core radio luminosity at 5 GHz. The 11 BCGs that have multi-frequency radio 
data available on them are indicated with the S - F symbols, for steep (a > 0.4) and flat/inverted (—0.2 < a < 0.4) spectrum sources 
respectively; (4) Intrinsic rest-frame X-ray luminosity in the 2 — 10 keV band of the nucleus (note that all are upper limits, except for 
Cygnus A and M87); (5) 2MASS K band absolute magnitude of the bulge; (6) A/bh estimated from the K band bulge luminosity using 
Eq. [2] (7) Mbh if the object sat on the fundamental plane, applying an average correction factor of log AA/bh = 0.8 ±0.6 to the K band 
derived BH mass; (8) Eddington ratio of the nuclear X-luminosity, using the BH mass in Column 6; (9) Eddington ratio of the nuclear 
X-luminosity, using the BH mass in Column 7. (10) Notes concerning the radio luminosities: i) derived from VLA observations at 4.9 
GHz; ii) derived from ATCA observations at 5.5 GHz; iii) extrapolated from the core 1.4 GHz flux and core spectral index (a = 0.56; 
iGiacintucci et al1l201ll) ; iv) no detection at 4.9 GHz; v) private communications (A. Edge); vi) extrapolated from the 1.4 GHz, 8.46 GHz 
and 28.5 GHz (BIMA) flux densities; vii) determined from EMSS; viii) from MH2007. All errors are 2<r. a K band absolute magnitude 
derived from UKIRT. 
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lOgA Mr.1t 
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A1835 


0.2532 


41.07±0.02 s 


<42 


25±0.19 


-27.36±0.28 


9.5±0.4 


10.3±0.7 


-5.4±0.4 


-6.2±0. 


7 


(i) 


A2204 


0.1522 


40.81±0.02 


<42 


21±0.13 


-26.57±0.27 a 


9.2±0.3 


10.0±0.6 


-5.1±0.3 


-5.9±0. 


7 


(i) 


A1664 


0.1283 


40.51±0.02 F 


<41 


18±0.48 


-26.15±0.22 


9.1±0.3 


9.9±0.7 


-6.0±0.5 


-6.8±0. 


8 


(ii) 


RXCJ1504.1-0248 


0.2153 


41.15±0.05 F 


<42 


43±0.28 


-26.59±0.36 


9.2±0.3 


10.0±0.7 


-4.9±0.4 


-5.7±0. 


7 


(iii) 


RXJ0439.0+0715 


0.2300 


<39.18 


<4l 


76±0.51 


-26.31±0.36 


9.1±0.3 


9.9±0.7 


-5.5±0.6 


-6.3±0. 


8 


(iv) 


A2390 


0.2280 


41.95±0.02 


<42 


03±0.20 


-27.05±0.34 


9.4±0.3 


10.2±0.7 


-5.5±0.4 


-6.3±0. 


7 


(v) 


A0478 


0.0881 


39.85±0.03 s 


<41 


53±0.20 


-26.72±0.14 


9.3±0.3 


10.1±0.7 


-5.9±0.3 


-6.7±0. 


7 


(vi) 


PKS0745-19 


0.1028 


40.13±0.02 s 


<4l 


,54±0.55 


-26.87±0.18 


9.4±0.3 


10.1±0.7 


-5.9±0.6 


-6.7±0. 


9 


(v) 


A2261 


0.2240 


39.70±0.07 s 


<4l 


,58±0.73 


-27.35±0.20 


9.5±0.3 


10.3±0.7 


-6.0±0.8 


-6.8±1. 





(i) 


Z2701 


0.2151 


40.40±0.02 F 


<41 


61±0.47 


-26.26±0.34 


9.1±0.3 


9.9±0.7 


-5.6±0.5 


-6.4±0. 


8 


(1) 


RXJ1720.1+2638 


0.1640 


40.00±0.04 F 


<41 


90±0.35 


-26.71±0.22 


9.3±0.3 


10.1±0.7 


-5.5±0.4 


-6.3±0. 


8 


(i) 


RXJ2129. 6+0005 


0.2350 


40.81±0.02 s 


<42 


23±0.27 


-26.73±0.30 


9.3±0.3 


10.1±0.7 


-5.2±0.4 


-6.0±0. 


7 


d) 


Z3146 


0.2906 


40.02±0.03 s 


<42 


85±0.13 


-26.45±0.56 


9.2±0.3 


10.0±0.7 


-4.4±0.3 


-5.2±0. 


7 


(i) 


MS1455.0+2232 


0.2578 


40.34±0.05 F 


<42 


29±0.19 


-27.13±0.28 


9.4±0.3 


10.2±0.7 


-5.3±0.4 


-6.1±0. 


7 


(i) 


MS2137.3-2353 


0.3130 


40.20±0.02 


<42 


,88±0.26 


-26.77±0.30 


9.3±0.3 


10.1±0.7 


-4.5±0.4 


-5.3±0. 


7 


(vii) 


Centaurus 


0.0104 


39.10±0.02 


<39 


.41±0.12 


-26.14±0.04 


9.1±0.2 


9.9±0.6 


-7.8±0.3 


-8.6±0. 


7 


(viii) 


Cygnus A 


0.0561 


41.43±0.02 


44.36±0.02 


-26.73±0.12 


9.3±0.3 


10.1±0.7 


-3.1±0.3 


-3.9±0. 


7 


(viii) 


M87 


0.00436 


38.88±0.02 


40.55±0.04 


-25.55±0.04 


8.8±0.2 


9.6±0.6 


-6.4±0.2 


-7.2±0. 


6 


(viii) 



of KFC2006, GCM2009 and PMK2012 are also included in 
our analysis of Section 3, and we show that our results are 
still consistent when using these correlations. 

2.2 Brightest Cluster Galaxies 

Our i nitial sample of BCGs includ e the 19 clusters stud- 
ied in iHlavacek-Larrondo fc Fabianl (|201ll ) and consists of 
18 clusters that have no detectable central X-ray nucleus, 
as well as Cygnus A which has clear non-thermal nuclear 
emission in the for m of an absorbed power law. We use the 
values obtained in IHlavacek-Larrondo fc Fabianl (|201ll ) for 
Z/x, all of which were derived from Chandra X-ray observa- 
tions. To obtain the values, we converted background sub- 
tracted number of counts of the nuclear region, within a 1" 
radius equivalent to the Chandra point s pread functio n, into 
fluxes using the web interface PIMM93 (|Mukaill 19931 ). Here, 
we include Galactic absorption and model the emission in 
the form of a power law with a spectral index equal to 1.9. 
Since the counts are dominated by thermal emission, with 
no clear evidence of non-thermal emission in the spectra, 
the fluxes we obtain are considered to be upper limits to the 
non-thermal contribution of the nucleus. The corresponding 
luminosities, along with the 2a errors are shown in Column 
4 of Tabled] 

The 5 GHz radio luminosities are obtained using a 
1 http:/ /heasarc. gsfc.nasa.gov/Tools/w3pimms. html 



combination of published values and archival maps, all of 
which have beam sizes < 5" (high angular resolution is 
needed to isolate the core emission). For RXC J1504. 1-0248, 
£sghz was derived using the 1.4 GHz core flux measure- 
ment (S'i.4GHz = 42.2 mjy ) and the core sp e ctral index 
(q = 0.56; S„ oc v~ a ) of IGiacintucci et all |201ll ). For 
A0478, we extrapolated the 5 GHz flux based on the peak 
Very Large Array (VLA) 1.4 GHz and 8.46 GHz fluxes (A 
configuration, project code AE117), as well as the 28.5 GHz 
measurement from the Berkele y-Illinois-Maryland Associa- 
tion (BIMA; ICoble et al.ll2007t ). The 5 GHz flux from the 
Einstein Observatory Ext e nded Medium-Sen sitivity Survey 
(EMSS; iGioia et al.lll990l ; IStocke et al.lll99ll ) was used for 
MS2137.3-2353. Finally, we adopted the 5 GHz core radio lu- 
minosities from lMerloni fc Heinzl (|2007l . hereafter MH2007) 
for the Centaurus cluster and Cygnus A. For the remaining 
objects, we computed the 5 GHz radio luminosities from the 
4.9 GHz C configuration VLA observations (project codes 
AE099, AE107, AE125; PI Edge) or 5.5 GHz Australia Tele- 
scope Compact Array (ATCA) observations (project code 
C1958), and use a conservative estimate by taking the peak 
intensity as the 5 GHz flux measurement. Our results are 
shown in Column 3 of Table [T] Errors are derived as the 
quadratic sum of the rms noise level in the radio map and 
the systematic error associated with the value. Systematic 
erro rs vary with freque ncy, but are on the order of 5 per cent 
(see ICarilli et al.l Il99ll ) . For simplicity, we therefore choose 
to compute the total error assuming a 5 per cent system- 
atic error and a 2a rms noise level. For the values taken from 
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20 22 24 26 28 30 32 34 
0.60log L x + 0.78log M BH 



Figure 1. Fundamental Plane of BH activity. We adopt the relation derived by MHD2003 and show the best-fitting regression with the 
solid line (see Eq.OJ. The sample of MHD2003 contains Galactic BHs (GBH), LINERs (L), transition nuclei (LINER/H II; T), Seyfert 
galaxies (Sy) and quasars (Q). We include the 18 BCGs from Table [l] (shown with the filled red and yellow circles), 16 of which have no 
detectable X-ray nucleus. For all of these objects, the BH masses were derived using the Mbh ~ Mk relation (Eq.[2|. We also include the 9 
BCGs from MH2007, 7 of which have BH masses derived from the Mbh — a relation (shown with the 5 pointed red and yellow stars). The 
remaining two (M87 and Cygnus A) have BH masses derived form dynamical modelling and are shown with the pink stars. Our results 
show that BCGs lie on or above the relation, i.e. are offset from the plane such that they appear too bright in the radio compared to their 
X-ray luminosities and predicted BH masses. A possible explanation is that the Mbh — °" ancl Mbh — Mk correlations systematically 
underestimate the masses of these BHs, implying that if these objects follow the FP, some would be ultramassive (Mbh > 1O 1O M0). 



the literature where no error estimate was available, we only 
considered systematic uncertainties and assume a 5 per cent 
uncertainty on the value. 

We then calculate the predicted BH masses for our 
B CGs using the M bh — Mk relation recently revised in detail 
bv lGrahaml (|2007D and shown in Eq.O 



log(M BH /M ) = -0.37 ±0 04 (M X + 24) + 8.29 



±0.08 



(2) 



We use K band Two Micron All Sky Survey (2MASS) 
magnitudes from the All-Sky Extend ed Source Ca t alogu e 
determined from fit extrapolation Jjarrett et al.l l2000t). 
which we correct for Galactic extinction ( Schlegel et al.l 
1 19981 ). We also correct for re dshift (Jf-correction) follo wing 
the latest IDL script^] from iBlanton fc Roweisl (|2007l ) and 
for evolution (_B-correction) using the K band 2MASS cor- 



http: / /howdy.physics. nyu.edu/index.php/Kcorrect 



rection estimated in lBell et al.l (|2003l . 0.8 x z). The resulting 
values and 2a errors are shown in Column 5 of Table [T] 

Three clusters did not have 2MASS observations 
(A2204, MACS J1532. 8+3021 and RX J1347.5-1145), but 
A2204 has a published United Kingdom I nfra- Red Te l escop e 
(UKIRT) absolute K band magnitude (|Stott et all 120081 ), 
where the authors estimate that the dominant error in their 
photometry is due to the fitting algorithm which can under- 
estimate the integrated brightness by up to 10 per cent of a 
magnitude. The error we show in Table[T]for A2204 is there- 
fore taken as 10 per cent of the value. Since the remaining 
two clusters (MACS J1532. 8+3021 and RX J 1347.5-1145) 
do not have 2MASS/UKIRT magnitudes, we do not include 
them in our final sample. 

The BH masses are then estimated using the absolute 
magnitudes (Mk) and Eq. [2] We use a Monte Carlo tech- 
nique to calculate the errors in the derived masses, where 
we assume that both parameters in Eq. [2] (0.37+0.04 and 
8.29+0.08), as well as Mx±AMj( are independent from one 
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Figure 2. Fundamental Plane of BH activity. Same as Fig. [TJ but 
focussed on the BCGs. Also highlighted are the 11 BCGs that 
have multi-frequency radio data available on them with the steep 
spectrum sources (a > 0.4) shown in black and the flat /inverted 
spectrum sources (—0.2 < a < 0.4) shown in light blue. 

another. We also assume that each follow a Gaussian distri- 
bution based on their values and associated errors. For each 
object, we proceed by selecting 100 random variables for all 
three distributions (0.37±0.04, 8.29±0.08 and M K ±AM K ), 
and run through all possibilities, each time computing the 
predicted BH mass. The final BH masses are then taken as 
the median values, and the 2a errors are calculated within 
the 2.2 and 97.6 percentiles (see Column 6 of Table [T]). 

Finally, we plot the location of each BCG on the FP 
following Eq. [Q in Fig. [j] Fig. [TJ also includes the 9 BCGs 
from MH2007: Cygnus A, NGC 1275 in Perseus, NGC 4486 
in M87, NGC 4696 in Centaurus, NGC 6166 in A2199, IC 
4374 in A3581, UGC 9799 in A2052, 3C 218 in Hydra A 
and 3C 388. In MH2007, the BH masses are derived from 
the Mbh — o relation or from dynamical mass measurements 
for M87 and Cygnus A. Note that Table [T] already includes 
NGC 4696 (Centaurus) and Cygnus A, and shows the BH 
masses predicted from the Mbh — Mk relation. We also in- 
clude M87 in Table [1] and derive its BH mass as predicted 
from the Mbh — Mk relation (from its 2MASS K band mag- 
nitude), in order to compare the predicted mass to the one 
measured from modelling of the kinematics. The X-ray and 
radio luminosities for M87 were taken from MH2007, where 
we assume a 5 per cent uncertainty on the first (equivalent 
to the systematic uncertainty in radio) and a 10 per cent 
uncertainty on the second. Finally, we include all the data 
points from MHD2003 in Fig.[TJ but remove Cygnus A, NGC 
4486 (M87), NGC 1275 (Perseus) and NGC 6166 (A2199), 
since they are already included in MH2007. 

In Fig. [21 we only show the location of the BCGs in the 
FP to illustrate more clearly where BCGs lie on the FP. Note 
that, most of our BCGs have no detectable X-ray nucleus 
(shown with the arrows), implying that they could lie even 
more leftward than shown in Fig. [2] 
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RX J0439. 0+0715 which has no detectable X-ray or radio 
nucleus. The BCGs are offset from the plane such that they 
appear too bright in the radio compared to their X-ray lu- 
minosities and predicted BH masses. Most of our objects 
lie in massive, strong cool core clusters that have power- 
ful outflows (I/ mcc h > 10 44-45 erg s _1 ) and no detectable 
X-ray nucleus ( Lx ^ 10 42 ergs -1 ). The majority of our 
objects are highly sub-Eddingtion with logA ;$ —5 (except 
for Cygnus A, see Table [T]). We therefore consider them to 
be in a state equivalent to the LH state of X-ray binaries. 

Although the scatter of the FP is large, it is unlikely 
that all BCGs lie randomly above the best-fitting relation; 
they should on average lie above and below the plane (i.e. 
not systematically above) . We therefore examine four differ- 
ent possibilities as to why BCGs appear to be systematically 
offset from the plane. 

3.1 Overestimating radio luminosities 

First, we investigate the possibility that we could be over- 
estimating the radio luminosities of all our BCGs. 

It is not clear why this would only be the case for BCGs 
considering that other AGNs such as Seyferts and quasars 
are also part of the MHD2003 sample and do not seem to 
be affected by such a bias. The 5 GHz radio luminosities 
of the MHD2003 sample were obtained from the literature. 
Their sample consists both of steep (a > 0.4) and flat (a < 
0.4) spectrum sources, and most of the luminosities were 
obtained from arcsec resolution VLA measurements and by 
integrating the fluxes. We have used a conservative approach 
to our 5 GHz luminosities, based on the peak fluxes and not 
integrated fluxes (peak fluxes are smaller than integrated 
ones). We should therefore not be overestimating the radio 
luminosities compared to the objects in MHD2003, yet our 
objects seem to lie above the FP. 

PMK2012 showed that the FP slopes change when steep 
spectrum sources are included. They also derive what they 
consider to be the most accurate FP regression to date (with 
the least scatter), based only on sub-Eddington accreting 
BHs with flat /inverted radio spectra, thus suggesting that 
only these sources should be included in the FP. Eleven of 
our objects in Table[T]have mult i- frequency radio data avail- 
able on them (mostly from the VLA archive), 6 of which have 
steep spectra (a > 0.4) and are shown in Column 3 of Table 
[T]with the s symbol. The remaining 5 have flat spectra with 
a varying between —0.2 to +0.4 and are shown in Column 
3 of Table [1] with the F symbol. Both steep and flat spec- 
trum sources are also shown in Fig. [2] with the black and 
blue points, respectively. If the offset observed for BCGs in 
the FP were due to the wrongfully included steep spectrum 
sources, we would have expected these sources to systemati- 
cally have the largest offsets. However, Fig.[2]shows that this 
is not the case, and that both steep and flat spectrum seem 
to spread out randomly in the scatter. It therefore seems 
unlikely that PMK2012's interpretation would explain why 
our objects are systematically offset from the plane. 



3 RESULTS 

Figs. [1] and [2] show that our BCGs systematically sit 
on or above the FP of BH activity. The only outlier is 



3.2 Underestimating X-ray luminosities: 
synchrotron cooling 

FKM2004 and KFC2006 argued that the X-ray emission 
for the most massive BHs (^1O 8 M ) could be strongly af- 
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fected by synchrotron cooling. Synchroton cooling is anti- 
correlated with BH mass, and for the most massive BHs, the 
cooling break occurs below X-ray wavelengths. Synchrotron 
cooling alters the spectral distribution above the cooling 
break and causes the X-ray luminosities to be underesti- 
mated from the true value for the most massive BHs. This 
makes BHs appear underluminous in the X-rays compared 
to their radio luminosities. We expect our objects to host at 
least 1O 8_9 M0 BHs. Synchrotron cooling could therefore ex- 
plain why our objects appear to be so radiatively inefficient 
and why they appear to be offset from the plane. 

PMK20f2 have explored this issue in more detail, and 
determined that even optical emission can be affected by 
synchrotron cooling in the most massive objects such as FR 1 
galaxies. They argue that spectral energy distribution (SED) 
modelling is necessary to obtain the luminosities and place 
the objects on the FP. Although this could explain why our 
objects are offset from the plane, SED modelling is beyond 
the scope of this paper, and we examine other possibilities 
as to why BCGs appear to lie above the FP of BH activity. 

3.3 Exceptional BHs in BCGs 

The third possibility is that BCGs are special in the sense 
that they occupy a particular place in the FP. This would 
imply that the BHs in these systems are operating in a differ- 
ent way than other BHs, which in itself is interesting. BCGs 
lie at the very centres of galaxy clusters and are surrounded 
by a substantial amount of hot dense gas. The most notice- 
able difference between BCGs and systems such as quasars 
and Seyferts is therefore the environment. It could be that 
the dense hot gas provides a unique environment that makes 
BCGs intrinsically more radio luminous for a given X-ray 
luminosity and BH mass, if the radio emission is due to syn- 
chrotron emission. 

3.4 Underestimating BH masses from the Mbh o 
and Mbh — Mr relations 

Finally, we examine the possibility that the Mbh — cr and 
Mbh — Mk correlations systematically underestimate the 
BH masses in BCGs. We base this idea on the recent studies 
that have obtained direct BH masses from dynamical mod- 
elling of the kinematics for a handful of BCGs (Macchetto 
et al. 1997; Houghton et al. 2006; Dalla Bonta et al. 2009; 
McConnell et al. 2011a,b), many of which seem to indicate 
that direct BH mass measurements from dynamical mod- 
elling are larger than the values predicted from the Mbh — cr 
and Mbh - M K correlations (e.g. M87, A1836-BCG, NGC 
3842 and NGC 4889). 

We do not expect all BCGs to intersect exactly the 
plane. However, it is unlikely that they randomly all lie 
above the plane in Fig. [2] Instead, they should on average 
lie above and below the plane. 

To illustrated this, we determine the average mass offset 
needed (log AA/bh) for BCGs to be consistent with the FP 
such that on average, they satisfy Eq. [3] log AA/bh therefore 
represents the offset needed so that BCGs lie on average 
above and below the plane. 

logisGH, = 0.601ogL x + 0.781og(MBH K ,„ + AM BH ) + 7.33 (3) 
All the clusters in Table [T] have error measurements for 



Table 2. Average BH mass offset from the FP - (1) Offset for 
the 18 BCGs that have a mass estimate from the Mbh — Mk 
correlation (Table [TJ. The errors are 2a and were determined 
using a Monte Carlo technique, see Section 3.4. (2) Probability 
that the mass offset is larger than zero in terms of per cent (and 
a). (3) Rough calculation of the average offset including the 18 
BCGs in Table Q] as well as the 7 BCGs from MH2007 that have 
a mass estimate from the Mbh — u correlation. 





(1) 


(2) 


(3) 




log AMjh,k 


Prob(> 0) 


log AM B h,k-<t 


MHD2003 


0.8±0.6 


> 99.6 (£3<t) 


1.0 


KFC2006 [MHD] 


1.5±0.7 


> 99.9 (> 3cr) 


1.6 


[FKM] 


2.1±0.7 


> 99.9 (> 3cr) 


2.3 


GCM2009 


0.5±0.6 


> 94.6 (>2<r) 


0.7 


PMK2012 


2.6±0.7 


> 99.9 (> 3cr) 


2.8 



Lsghz, Lx and Mbh, A" (from K band magnitudes), thus 
allowing us to determine the average mass offset needed by 
using a Monte Carlo technique. Here, we assume that Lsghz, 
Lx and Mbu,k are independent, and that each follow a 
Gaussian distribution based on their values and associated 
uncertainties. For each of the 18 BCGs in Table [1] we as- 
sign 500 random variables to the Lsghz distribution, 500 
to the Lx distribution and 500 to the ALbh, k distribution. 
Using 1000 variables instead of 500 for each distribution 
does not change our results. Then, for each of the 500 3 pos- 
sibilities and for each of the 18 objects, we calculate the 
log AMbh, k needed for the object to satisfy Eq. [3] For each 
of the 500 3 possibilities, we then calculate the average value 
of log AA/bh, k over the 18 objects. This allows us to build a 
distribution containing 500 3 estimates of the average offset 
needed for BCGs to be consistent with the FP. Our results 
do not change significantly if median values are used instead 
of average values. 

The final mass offset is taken as the median value in the 
log AA/bh, k distribution and the 2a errors are taken within 
the 2.2 and 97.6 percentiles. We repeat this calculation for 5 
different FP regressions: MHD2003, both for the revised re- 
lations of FKM2004 and MHD2003 in KFC2006, GCM2009 
and PMK2012. The results are shown in Column 1 of Table 
[2] We emphasize that 16 of the 18 BCGs in Table Q] have no 
detectable X-ray nucleus. If we were to remove the 2 BCGs 
with X-ray nuclei from the calculations, we would obtain 
the same mass offsets. The calculations are therefore domi- 
nated by the non-detections, and the mass offsets should be 
regarded as the minimum offset needed for BCGs to lie on 
the FP. 

We have also considered the 7 BCGs from MH2007, 
which have BH mass measurements based on the Mbh — <r 
relation. Since MH2007 do not include uncertainty measure- 
ments on their values, we only compute the average mass 
offset needed for all BCGs to be consistent with the FP 
such that they satisfy Eq. [3] (no Monte Carlo calculations 
performed over the uncertainties). Here, all 18 BCGs from 
Table Q] and the 7 BCGs from MH2007 are included in the 
calculation. Note that NGC 4696 (Centaurus) is counted 
twice in the calculation, once for the mass derived from the 
K band magnitude and once for the mass derived from the 
velocity dispersion (a). The offsets we find are shown in Col- 
umn 3 of Table [2] and represent the offsets needed such that 
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35 40 45 50 55 

log L x + 0.88log M BH 



Figure 3. Same as Fig.[T] but where we consider the best-fitting relation of PMK2012 (black line) which is considered to be the most 
robust to date. It is also the relation that has the largest mass offset for BCGs from the FP (see Tabled . Shown are the data points used 
in PMK2012 to derive the best-fitting relation, which only include sub-Eddington BHs with flat/inverted radio spectra (GBHs, Sgr A*, 
LLAGN, as well as HBLs from SDSS). The BCGs are illustrated with the same symbols as in Fig. [2]and are offset from the plane such 
that they appear to be less massive than predicted from their X-ray and radio luminosities. A possible explanation is that the Mbh — cr 
and Mbh ~~ Mjc correlations systematically underestimate the masses of BHs in BCGs. 



all 25 BCGs lie on average above and below the best-fitting 
plane. 

Table[2]shows that all of the offsets are positive, i.e. con- 
sistent with the idea that the A/bh — a and A/bh — Mjc cor- 
relations systematically underestimate BH masses in BCGs. 
Interestingly, we find that the PMK2012 regression, consid- 
ered to be the most robust to date, has the largest mass 
offset. Since this is currently the most accurate FP regres- 
sion, we show the location of our BCGs in this plane in Fig. 
E2 Note that, PMK2012 only considered sub-Eddington BHs 
(i.e. BHs considered to be in the LH state) with flat/inverted 
radio spectra. By limiting the sample to these sources, 
they obtained the most tightly constrained regression co- 
efficients to date. Their sample consists of GBHs, Sgr A*, 
LLAGN and 39 BL Lac objects from the Sloan Digital Sky 
Survey (SDSS). The latter are relativistically beamed sub- 
Eddington AGN where the jet core emission can be better 
isolated due to the beaming. The sample of BL Lac ob- 
jects only includes those that have peak synchrotron lumi- 
nosities in the soft X-rays, referred to as high-energy cut- 
off BL Lac objects (HBL). In this case, the X-ray emission 



is mostly considered to be synchrotron emission, whereas 
with low-energy cutoff BL Lac objects (LBL), the peak syn- 
chrotron luminosities occur at near-infrared wavelengths, 
and the synchrotron emission is most likely contaminated 
with synchrotron self-Compton and external inverse Comp- 
ton emission. Since the FP is thought to arise from BHs in 
the LH state dominated by synchrotron jet emission, it is 
necessary to only consider the sample of HBL as opposed 
to LBL. In Fig. g]we show PMK2012's best-fitting relation 
along with their sample of 82 objects. This figure shows 
clearly that our BCGs lie significantly offset from the plane, 
such that they appear to be less massive than that predicted 
from their X-ray and radio luminosities. 

Table [2] also shows that the best-fitting relation by 
GCM2009 is only consistent with a positive offset to a ~ 95 
per cent confidence level. GCM2009 only uses objects which 
have direct BH mass measurements from dynamical mod- 
elling. Additionally, they were able to determine a best- 
fitting relation that estimates a BH mass from its X-ray and 
radio luminosities. If we use this relation, we find that the 
masses predicted by GCM2009 are on average larger than 
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those predicted by the Mbh — a and Mbh — Mk relations, 
although the scatter remains large. However, if we use the 
best-fitting relation they find when excluding Seyfert galax- 
ies, we find that the masses predicted by GCM2009 are sys- 
tematically larger than those predicted by the Mbh — cr and 
Mbh — Mk relations, for all 25 of our BCGs except for the 
BCG in A2261. Seyfert galaxies are thought to have higher- 
accretion rates, and are therefore not in a state equivalent 
to the low-accreting LH state of X-ray binaries. 



4 THE EXISTENCE OF UMBHS IN BCGS 

If BCGs truly follow the FP of BH activity and that our 
X-ray and radio core luminosities have not been underesti- 
mated or overestimated respectively, then our results imply 
that the intrinsic masses of the BHs are higher than those 
predicted from the Mbh — o and Mbh — Mk correlations. 

The sample of MH2007 contains 9 BCGs, 2 of which 
have precise BH mass measurements from dynamical mod- 
elling. The first is Cygnus A, which has a measured BH 
mass of 2.5 x 1O 9 M and a predicted BH mass from the K 
band relation of 2 x 1O 9 M0, both of which are consistent 
with one another. However, the measured BH mass for M87 
is significantly higher (3 x 10 9 M Q , MH2007; 6.6 x 10 9 M o , 
iGebhardt et al.lboill ') than that predicted from the K band 
relation (0.8 x 10 9 Mq). The offset is such that it agrees with 
our calculations. 

For the remaining 7 BCGs in MH2007, our results imply 
that the true BH masses lie between 1O 9 M and 10 1C, M Q . 
However, for the 16 BCGs in Table [1] that have no de- 
tectable X-ray nucleus and that lie in some of the most ex- 
treme clusters of galaxies (i.e. some of the most massive and 
strong cool core), if we apply the average mass increase of 
logAA^BH = 0.8, then our results imply that the true BH 
masses lie between 8 x 10 9 M Q and 2 x 10 10 M Q . Further- 
more, if we apply the average mass increase derived when 
including the BCGs from MH2007 (logAM BH = 1.0, see Ta- 
ble 2), then the majority of these BHs are ultramassive with 
M BH = (1 -4) x 10 10 M o . 

The existence of UMBHs in BCGs has recently been 
confirmed by McConnell et al. (2011a) who find a 2.1 x 
10 10 M o BH in NGC 4889. Yet, NGC 4889 lies in a fairly av- 
erage cluster with no cool core associated with it. It therefore 
does not require an active nucleus to prevent the surround- 
ing gas from cooling. On the other hand, the central BHs 
in massive and strong cool core clusters with Lx > 10 45 
erg/s and t coo \ < 3 Gyr (such as the majority in Table [T]) 
must be injecting extreme mechanical energies into their 
surrounding medium to prevent the gas from cooling, on 
the order of 10 44 " 45 erg s" 1 . Out of all BCGs, these are 
where the BH must have accreted a substantial amount of 
mass to power the outflows (> 10 9 M3). Since these BCGs 
also lie in the most massive clusters, and BCG mass scales 
with cluster mass, they should host the most massive BHs 
compared to other BCGs. The most massive BHs should 
therefore reside in these massive and strong cool core clus- 
ters, and our results support this claim, while predicting 
that many will have UMBHs at their centres. If confirmed, 
our results will have important ramifications fo r the forma- 
tion a nd evolution of BHs across cosmic time (lLauer et al.l 
l2007i ; iMcNamara et all 120091 ; iNataraian fc Treisterl 120091 ; 



Hlavacek-Larrondo fc Fabianl 1201 ll ; IMcNamara et al.l [ioill ; 
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It is also important to note that BCGs have undergone 
the most dramatic environments, subject to major mergers 
in the past and extreme AGN outflows for the past several 
Gigayears. It is not surprising then, that they would have 
had the opportunity to grow to such masses. There are two 
possible scenarios in which BHs can grow to be ultramas- 
sive. The first is through hiera rchical mergers, as supported 
by numerical calculations fe.g. lYoo et al.|[2007l . who predict 
the existence of a rare population of UMBHs in the local 
Universe). The second is from high redshift "seeds" and is 
based on the observation tha t quasars exist fr om as early as 
z of about 6 |VikhliniiJl20il ; lFan~ et al .1120061 ). UMBHs can 
therefore form from these high redshift "seeds", and evolve 
into present-day BCGs, which are the most massive galax- 
ies of the local Universe. INataraian fe Treisterl |2009l ) argue 
however that although UMBHs may exist, the maximum 
mass they can reach is around ~ 10 10 Mq. 



5 CONCLUDING REMARKS 

We have identified two possibilities as to why BCGs seem 
to be systematically offset from the FP of BH activity, as- 
suming that our X-ray and radio core luminosities are cor- 
rect. The first is that the BHs in BCGs are supermassive 
(Mbh ~ 1O 9 M0), but operate differently from other BHs 
and lie in a particular place in the FP of BH activity. The 
second is that these BHs follow the FP of BH activity but 
not the standard Mbh — o and Mbh — Mk correlations, 
thus predicting that many of the BHs in massive and strong 
cool core clusters are ultramassive (Mbh > 10 10 Mq). Our 
results therefore carry significant implications for the forma- 
tion and evolution of BHs, as well as the connection between 
the central BH and its host galaxy. Only by obtaining direct 
BH masses from dynamical modelling of the kinematics, as 
opposed to relying on known correlations, can we determine 
if these extreme BHs are truly ultramassive. 
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